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ABSTRACT 

We present wide field near-infrared photometry of 12 Galactic globular clusters, typically extending 
from the tip of the cluster red giant branch (RGB) to the main sequence turnoff. Using recent 
homogenous values of cluster distance, reddening and metallicity, the resulting photometry is directly 
compared to the predictions of several recent libraries of stellar evolutionary models. Of the sets 
of models investigated, Dartmouth and Victoria-Regina models best reproduce the observed RGB 
morphology, albeit with offsets in J — Ks color which vary in their significance in light of all sources 
of observational uncertainty. Therefore, we also present newly recalibrated relations between near- 
IR photometric indices describing the upper RGB versus cluster iron abundance as well as global 
metallicity. The influence of enhancements in alpha elements and helium are analyzed, finding that 
the former affect the morphology of the upper RGB in accord with model predictions. Meanwhile, 
the empirical relations we derive are in good agreement with previous results, and minor discrepancies 
can likely be attributed to differences in the assumed cluster distances and reddenings. In addition, 
we present measurements of the horizontal branch (HB) and RGB bump magnitudes, finding a non- 
negligible dependence of the near-IR HB magnitude on cluster metallicity. Lastly, we discuss the 
influence of assumed cluster distances, reddenings and metallicities on our results, finding that our 
empirical relations are generally insensitive to these factors to within their uncertainties. 

Subject headings: globular clusters: general — globular clusters: individual(NGG 104, NGG 288, NGG 
362, NGG 1261, NGG 1851, NGG 2808, NGG 4833, NGG 5927, NGG 6304, NGG 
6496, NGG 6584, NGG 7099) — stars: infrared 


1. INTRODUCTION 

Galactic globular clusters (GGGs) play a critical role 
as templates of evolved stellar populations, functioning 
as testbeds for stellar evolutionary models. At opti¬ 
cal wavelengths, huge advances in our understanding of 
GGGs have been made via extensive photo metric surveys 
of large samples of GGGs both from space (IP lotto et al.l 


lar ge __ __ 

l20Q^ : ISaraiedini et al.ll20(y^lPiotto et al.l[2014fl. and from 

the ground (e.g. iRDsenberg et al.l 119991: iStetsonI l2QQQf ). 

The homogeneity of these databases has facilitated a va¬ 
riety of comparisons between prec ise photometry and 
existi ng evolutionary models (e.g. iMarin-Franch et al.l 
l2QQ9l : Dotter et al. 2 010, hereafter iDlOl : VandenBerg 
et al. 2013, hereafter IV13f ). However, at near-infrared 
(near-IR) wavelengths, such comparisons remain few de¬ 
spite a growing database of observations. iBrasseur et al.l 
(|2010D compared VJKs photometry of 6 GGGs (and 
the old open cluster NGG 6791) to the predictions of 
the l atest Victoria-Regina models (jVandenBerg et all 
120141) . finding that the models fail to reproduce the 
observed morphology of cluster RGBs, at least at 
lower metallicities. On the other hand, the Victoria- 
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Regina models seem to function well in the near-IR 
when compared to M4 and NGG 6723 (iHendricks et al 


2012 ) and the main sequence of NGG 3201 (jBono et al 


20101 ). Other direct comparisons between near- 
IR GGG isochro n es and p hotometry were performed 
bv IValenti et al.l (l2004al|bl) u s ing t h en-recent models 
(iCaloi. D’Antona fc Mazzitellil Il997l: l^laris fc Gassisil 


119971 : iStraniero. Ghieffi fc Limongil 19971) a nd bolometidc 

correc tions (jMontegriffo et al.l 1199^ . and iSalaris et al.l 
(j2007l ) compared the observed HB and RGB bump and 
lum inosity function of 47 Tn c to a-enhanced BaSTI mod¬ 
els (jP letrinferni et aP 120061 ). The lack of existing com¬ 
parisons between isochrones and data in the near-IR 
is not due to a lack of model predictions, and in fact 
ISalaris fc Girardil (j2002[ ) tabulate the predicted mean ab¬ 
solute K magnitude Mk of cluster HBs as a function of 
age and metallicity. However, at ages typical of GGGs 
(>10 Gyr) these predictions have never been tested thor¬ 
oughly as the combination of sufficiently accurate ages 
and iF-band photometry was lacking. 

The systematic observation of GGGs in the near-IR has 
been undertaken largely by E. Valenti and collaborators, 
yielding a database of near-IR photometry of optically 
well-studied GGGs, including relatio ns between observed 
photometric features a nd metallicity ([Ferraro et al.[l2QQQl: 
IValenti et al.l 2QQ4allbl: iFerraro. Valenti fc Origlial I2QQ6D 


m add ition to lCho fcLed ( 2QQ2D and llvanov fc Borissoval 
(|2QQ2f ). who undertook similar analyses using photom- 
etry from the Two Micron All Sky Survey (2MASS; 
iSkrutskie et al.l [20061 ). These studies have since served 
as a de facto template for old (>10 Gyr) stellar popu- 
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lations in the near-IR, and have therefore seen numer¬ 
ous applications. These include the measurement of 
photometric distances, reddenings and metallicities for 
many poorly stud i ed GGCs located towards the Galactic 
bulge (jKim et al.l I2QQ6I: iChnn et all l2QlQl : IValenti et al. 
2Q1C ), Local Group d warf galaxies (e.g. iGorski et al. 
2(nl [HdrG~ZI l2(^ and redd ening and metahicity 
maps of the Galactic bulge (e.g. iGonzalez et al.l l2QllL 

Imi). 

Meanwhile, the proliferation of space-based photomet¬ 
ric GGG surveys has yielded homogenous databa ses of 
GGC cluster distances, reddenings and ages (e.g. Eia; 
IVI while high resolution multi-object spectroscopic 
campaigns have produced detailed chemical abundance 
studies which include large homogenous GGG abundance 
databases of [Fe/H] (Garretta et al. 2 QQ9a, hereafter 

a as well as light a nd a-elements (e.g iGarretta et all 
[2oomi^ . [201Q313). With the goal of leveraging 
these improvements together, we present near-IR pho¬ 
tometry of a sample of optically well-studied GGGs, 
which we employ to reexamine relations between pho¬ 
tometric features and metahicity. We then compare our 
observations to the predictions of several sets of evolu¬ 
tionary models in near-infrared colors for the first time. 

In the next section, we describe the observations and 
data reduction. In section three, we present color- 
magnitude diagrams and empirical fiducial sequences for 
all of our target clusters. Section four contains a compar¬ 
ison between the observed RGB fiducial sequences and 
predictions of five sets of evolutionary models, as well as 
updated calibrations of photometric indices versus clus¬ 
ter metahicity. In section five, we present observed mag¬ 
nitudes of the horizontal branch and red giant branch 
bump, as well as empirical relations between the near- 
IR absolute magnitude of these features versus cluster 
metahicity. Section six consists of a discussion of current 
uncertainties in cluster distance, reddening and metahic¬ 
ity, and their inhuence on our empirical calibrations, and 
in the final section we summarize our results. 

2. OBSERVATIONAL DATA 
2.1. Observations and Pre-Processing 

Observations of our 12 target clusters were obtained 
with the Infrared Side Port Imager (ISPI) mounted on 
the 4m Blanco telescope at Gerro Tololo Inter-American 
Observatory. The HAWAII-2 2048x2048 pix detector 
has 0.305"/pixel, giving a field of view 10.25 arcmin per 
side. Imaging was obtained in the J and Ks filters over 
the course of three runs between 2008 and 2010, with me¬ 
dian seeing ranging between 0.8"and 1.4", and a log of 
the observations is given in Table [TJ In order to mitigate 
the effects of persistence, saturation, and cosmetic de¬ 
fects, a two-step (for the 2008 run) or five-step (ah other 
runs) dither pattern was employed for each cluster, where 
each individual 60s exposure in the dither pattern is com¬ 
prised of 3-12 coadds (as given in the last two columns 
of Table [T]) to optimize dynamic range. As near-infrared 
imaging is typically limited by the sky background, which 
can vary both spatially (on scales smaller than the detec¬ 
tor field of view) and temporally (on timescales of min¬ 
utes), careful subtraction of this sky background is crit¬ 
ical to mitigating photometric errors. Therefore, each 
on-target dither sequence for each cluster was interca¬ 


lated with identical sequences targeted at an offset sky 
field located about ^15' away. These offset sequences 
were used to create a combined sky frame corresponding 
to each on-target sequence. Preprocessing, including the 
construction of bad pixel masks, bias subtraction, and 
hat helding, was accomplished using IRAhl^ tasks cus¬ 
tomized for ISPI within the cirred packagqj. Next, sky 
subtraction was optimized by using a multi-step itera¬ 
tive procedure: Sources are rejected to ht an initial sky 
background, before median combining and subtracting 
these individual background frames from the original sky 
frames to generate more hne-scale sky frames. These im¬ 
proved sky frames are again median combined and added 
to the initial background, yielding a hnal sky frame for 
each exposure sequence. The high-order (4 <n<6) dis¬ 
tortion terms present in ISPI images were corrected us¬ 
ing the IRAF ccmap and ms c image tasks by matching 
many (typically several hundred) well-detected sources 
to 2MASS. Finally, a low-order background is ht to each 
sky-subtracted, distortion corrected image to remove any 
residual large-scale gradients. 

2.2. Photometry and Calibration 

Point-spread function htting (PSF) photom¬ 
etry was performed via ite rative usage of the 
DAOPHOT /ALLFRAME suite (fStetsoulIT9871 . 1199^ as 
described in iManro et al ] (l2nT^ . To limit our analysis 
to well-measured point sources, the resulting instru¬ 
mental catalogs were culled to retain only stars with a 
photometric error a{J — Ks) <0.1 and |sharp| <0.2. 

This latter cut was found to efficiently reject both 
spurious detections as well as stars with colors substan¬ 
tially ahected by crowding, illustrated in Fig. [T] for the 
example of NGG 362. There, detections which failed 
the sharpness cut, plotted in red, are found not only 
to cluster around bright sources spatially (as seen in 
the upper right panel), but also have RGB colors which 
are scattered preferentially blueward due to blending, 
seen in the color-magnitude diagram (GMD) as well 
as a color histogram of the lower RGB (bottom right 
panel of Fig. [T])» This effect is discussed, for example, by 
iBergbnsch fc StetsonI (|2QQ9l ). and demonstrates why a 
fairly stringent sharpness cut is needed despite the loss 
of some real sources in order to guard against systematic 
color offsets as a result of photometric blends. We 
caution that by electing straightforward self-consistent 
photometric quality cuts, a small fraction of spurious 
detections may remain in our final catalogs, although 
this has a negligible impact on our results since our 
analyses typically employ stars well above our detection 
limit, and any remaining spurious detections likely fall 
well redward of the cluster sequences (see Fig. [T]). 

In order to facilitate comparison and concatenation 
with existing GGG near-IR photometry (see Sect. Ira 
as well as direct comparis on t o multiple sets of evolu¬ 
tionary models (see Sect. lOTi . we have chosen to cal¬ 
ibrate our instrumental catalogs to the 2MASS photo¬ 
metric system. Stars astrometrically matched to 2MASS 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ See http : //www. ctio. noao. edu/noao/content/ISPI-Data-Reduction 
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TABLE 1 

Log of ISPI Observations 


Date 

Gluster 

N(J) 

N(Ks) 

Goadds( J) 

Coadds {Ks) 

13 Aug 2008 

NGG 104 

4 

12 

12 

6 


NGG 6496 

6 

29 

6 

12,6 

30 Sep 2009 

NGG 1851 

4 

29 

3 

6 


NGG 288 

10 

26 

3 

6 

01 Oct 2009 

NGG 362 

15 

34 

3,6 

6 


NGG 1261 

10 

29 

6 

6 


NGG 7099 

15 

34 

3,6 

6 

28 Apr 2010 

NGG 2808 

10 

34 

4 

6 


NGG 6304 

10 

35 

4 

6 

29 Apr 2010 

NGG 4833 

5 

40 

4 

6 

30 Apr 2010 

NGG 5927 

10 

33 

4 

6 


NGG 6584 

9 

34 

4 

6 



- 0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.3 0.4 0.5 0.6 0.7 

(J-Ks) (J-Ks) 


Fig. 1.— Left: CMD of NGC 362 showing all stars which passed 
our photometric quality cuts (black) as well as those which failed 
the sharpness cut (red). The dashed horizontal lines indicate the 
magnitude range used to construct a color histogram of the lower 
red giant branch (lower right), illustrating the bias caused by stars 
with higher absolute values of the sharpness parameter. The spa¬ 
tial distribution of the two samples relative to the cluster center 
(upper right) illustrates the efficiency of the sharpness criterion to 
eliminate spurious detections. 


were only used as photometric calibrators if they sat¬ 
isfied a more stringent set of criteria: They must be 
matched to a 2MASS point source to within the astro- 
metric rms (<0.2"), isolated (lacking neighbors within 4 
mag inside a 2.5" radius, corresponding to a contami¬ 
nating flux of <2.5%) and be bright enough to remain 
unaffected by crowding in 2MASS (see below). To calcu¬ 
late transformations from the instrumental magnitudes 
in our PSF catalogs to the 2MASS photometric system, 
we employed classical linear transformation equations of 
the form m — M = a + b{J — Ks), where m and M 
denote instrumental and standard magnitudes respec¬ 
tively. We solve for the coefficients a and b using least 
squares fitting, but employing weighting factors to down¬ 
weight discrepant data points in lieu of a sigma clip0 


TABLE 2 

Photometric Calibration Uncertainties with Respect to 
2MASS 


Gluster 

wrms 

a(zpt) 

N 


J 

Ks 

J 

Ks 


NGG 104 

0.027 

0.034 

0.005 

0.006 

647 

NGG 288 

0.019 

0.024 

0.013 

0.017 

89 

NGG 362 

0.025 

0.031 

0.007 

0.010 

222 

NGG 1261 

0.035 

0.034 

0.014 

0.014 

76 

NGG 1851 

0.026 

0.027 

0.010 

0.011 

118 

NGG 2808 

0.024 

0.029 

0.006 

0.008 

347 

NGG 4833 

0.023 

0.026 

0.009 

0.010 

164 

NGG 5927 

0.028 

0.032 

0.005 

0.006 

495 

NGG 6304 

0.032 

0.033 

0.010 

0.011 

360 

NGG 6496 

0.022 

0.029 

0.006 

0.007 

220 

NGG 6584 

0.025 

0.030 

0.007 

0.009 

135 

NGG 7099 

0.026 

0.025 

0.033 

0.036 

61 


(|Manro et al.l l2Ql^ iCohen et al.l l2Ql4 ). and we apply 
these same weighting factors to calculate the weighted 
root mean square deviation (wrms) of the residuals. A 
comparison between the results of our calibration proce¬ 
dure and the magnitudes reported in the 2MASS PSC are 
shown in Fig. El with the wrms given in each plot. The 
asymmetry seen at the faint end among stars positionally 
matched to 2MASS is due to the brighter completeness 
limit and lower spatial resolution of 2MASS, resulting 
in a systematic brightward deviation in 2MASS magni¬ 
tudes. This effect serves as a cautionary note against 
calibrating to 2MASS in cases where only the faintest 
2MASS stars are unsaturated, and for this reason we 
exclude faint 2MASS sources {Ks ^14 depending on 
stellar density) from use as photometric calibrators (see 
Fig. Ej)- Statistics regarding calibration uncertainties for 
each cluster, including the wrms values, photometric zero 
point uncertainties, and the number of stars from 2MASS 
used for photometric calibration, are summarized in Ta¬ 
ble El As our photometry saturates slightly below the 
tip of the cluster red giant branch in some cases, we have 
supplemented our catalogs with 2MASS photometry for 
bright stars which are absent due to saturation. 

2.3. Comparison to Previous Photometry 

Four of ou r targe t clust ers have been observed by 
IValenti et al.l (j2QQ4ai l2QQ5[ lFI. and near-IR photometry 


^ The algorithm is based on a series of five lectures pre- http://www.cade.hia.nrc.gc.ca/community/STETSON/\homogenous/Technique 

sented at ”V Escola Avancada de Astrofisica” by Peter Stetson ^ Also see the Bulge Globular Gluster Archive at 

http://ned.ipac.caltech.edu/level5/Stetson/Stetson_contents. htmhttp ://www.bo.astro.it/$\mathaccent"707E\relax{}$GC/ir_archive 
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Fig. 2.— A comparison of our calibrated photometry in the J (left) and Ks (right) filters with photometry from the 2MASS PSC. All 
matched stars are shown in grey, and 2MASS stars used for photometric calibration are overplotted in black. In each plot, the solid grey 
line represents equality while the dashed grey lines represent the weighted RMS deviation among the calibrators. 


of 47 Tuc has been presented by iSalaris et al.l (|2QQ7[ ). 
allowing a direct comparison of our 2MASS-calibrated 
JKs photometry to theirs. In all cases, we recover the 
majority of previously detected sources, and a cluster- 
by-cluster comparison of magnitude difference as a func¬ 
tion of magnitude is shown in Fig. [3l There, we have 
computed the mean magnitude difference for each clus¬ 
ter in each filter using a weighted 2.5(7 clip in magnitude 
bins while excluding stars faintward of the observed lu¬ 
minosity function peak. Our calibrated photometry is in 
agreement with earlier studies in light of total calibra¬ 
tion uncertainties, and although offsets of >0.05 mag are 
seen for NGC 288 in the Ks band (and to a lesser extent 
NGC 7099 in J), the direct comparison to 2MASS in 
Fig. Ogives no reason to be doubtful about the calibra¬ 
tion of these clusters given our calibration uncertainties 
listed in Tabl e [2] and the zeropoin t uncertainty of ±0.05 
mag cited by I Valenti et al.l (j2005f ) . 

In addition, there is one cluster in our sample (NGG 
1851) in common with the study of iBrasseur et al.l 
(j2010[ ). To compare their {V — Ks) fiducial sequence 
with our observations, we have matched our photom¬ 
etry to publicly available optical photometry from the 
archive of P. B. StetsorrM as well as the photomet¬ 
ric c atalog available from t he AGS GGG Treasury Sur¬ 
vey (|Sarajedini et al.ll2007D. In th e latter case, the off- 
sets listed by iHempcl et al] (j2014D were applied to the 
iSirianni et al.l (j^05f) transformed magnitucies to place 
them on the photometric system employed by Stetson. 
As the AGS catalogs are astr ometrically calibrated to 
2MASS (j Anderson et al.ll2QQ8f ). a simple matching algo- 
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Fig. 3. — C omparison of our 2M ASS-calibrated JKs photometry 
with that of ISalaris et alJ (120071) for N GC 104 (top panel) and 
I Valenti et al.l ii2nn4an and I Valenti et all (j 20051 1 for the remainder of 
the clusters illustrated. The solid grey line represents equality, and 
the dotted grey line represents the mean offset, calculated using a 
2.5(7 clip in bins of 1 mag and given in the upper left of each panel. 


rithm with an initial tolerance of 0.6" and an iterative 
5(7 clip functions well, and we recovered 99% of ISPI 
sources in the AGS field of view with an astrometric rms 
10 See http://www3. cade-ccda.hia-iha.nrc-cnrc.gc.ca/en/conrniuni^/^Tp^imiArhe resulting matched CMD is shown in Fig.H] 
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Fig. 4. — Observed Ks, (V — Ks) CMD of NGC 1851, obtained 
by matching our ISPl photometry with publicly available opti- 
cal photometry (see text for details). The fiducial sequence of 
IBrasseur et al.l (|2010l') is overplotted as a black solid line. 

with t he optical-IR fiducial sequence of IBrasseur et all 
(|2Q1Q[) over plotted. 

3. COLOR-MAGNITUDE DIAGRAMS AND FIDUCIAL 
SEQUENCES 

3.1. Construction of Fiducial Sequences 

To construct fiducial sequences representative of the 
RGB in each cluster, we fit a low-order (n <3) polyno¬ 
mial in the (J — Ks) vs. K s plane in a similar man¬ 
ner as previous st udies (e.g. iFerraro et al.l 1 1 9991 I2QQQI: 
iCohen et al.ll2Ql3) : First, a rough color-magnitude cut 
is made in the cluster CMD to select the region of the 
RGB. Next, the RGB is divided into bins of 0.5 mag, and 
the median color and magnitude (and their uncertainty) 
are calculated in each bin. A polynomial is then fit to 
these median points, iteratively rejecting stars more than 
2cr from the polynomial in each bin until convergence is 
indicated by the total number of stars remaining constant 
to within 2% since the previous iteration. 

In Fig. Owe present J — Ks^ Ks CMDs of all of our tar¬ 
get clusters, with the RGB fiducial sequences overplotted 
in red and median photometric errors in bins of 1 mag in¬ 
dicated along the right-hand side of each CMD. All vari- 
bles matched to our photom e try fr om the most recent 
version of the I Clement et al.l (|2Q1Q[ ) catalog of variable 
stars in GGC0 are overplotted as blue diamonds and 
excluded from further analyses. 

Our photometric catalogs are available electronically, 
and their format is illustrated in Table [3l 

3.2. The Absolute Magnitude Plane: Cluster Distance, 
Reddening and Metallicity 

One of the primary goals of the present study is to 
use a self-consistent set of cluster distances, reddenings 
and metallicities to compare observed evolutionary se¬ 
quences to models. Furthermore, with an eye towards 

See http :// www.astro.utoronto.ca/~ cclement/cat/listngc.html 


TABLE 3 

NGG 104 Photometric Gatalog 


RA (J2000) 

Dec (J2000) 

J 

(jJ 

Ks 

<jKs 

6.283727 

-72.069279 

10.9506 

0.0263 

10.3440 

0.0044 

6.283700 

-72.031622 

14.2151 

0.0155 

13.7408 

0.0075 

6.283753 

-72.133077 

15.2417 

0.0217 

14.9558 

0.0144 

6.283345 

-72.093668 

16.7742 

0.0621 

16.5855 

0.0298 

6.283288 

-72.071959 

17.8108 

0.0765 

17.2019 

0.0482 

6.283036 

-72.096528 

17.0105 

0.0611 

16.9107 

0.0429 

6.282974 

-72.112954 

16.2559 

0.0342 

16.1854 

0.0254 

6.282944 

-72.063674 

11.2960 

0.0234 

10.6157 

0.0042 

6.282786 

-72.060823 

16.6431 

0.0446 

16.2261 

0.0231 

6.282719 

-72.022186 

17.9398 

0.0723 

17.4222 

0.0490 


Note. — The full photometric catalogs for all target clusters are 
available electronically; a portion is shown here to illustrate their 
form and content 


heavily reddened stellar systems for which the use of op¬ 
tical photometry to measure cluster parameters is pro¬ 
hibitive, we also rederive empirical relations describing 
the RGB shape as a function of metalli city. To this end, 
we select distances and r edden ings from lDlOl and \F e/H] 
(and uncertainties) from IGQ9I . For consistency with IDIQI 
we empl oy the total t o selec tive extinction ratios Ry 
given by iSirianni et al.l (|2QQ5[) for a G2 star, and values 
of Aj/E{B — R)=Q.899 and AjcJE{B — 'F)=0.366 from 
iGasagrande fc VandenBer"g| (|2Q14D . For the two GGGs 
excluded from iDlQl due to the presence of multiple stel¬ 
lar populations (NGG 1851, NGG 2808), we employ pa¬ 
rameters obtained identically as for the remainder of the 
dToI sample (A. P otter, private communication, also see 
Milone et al.ll2Ql3) . 

Our adopted values for the cluster distance moduli, 
reddening and metallicity are listed in Table IH In light 
of cluster-to-cluster variations in a-enhancement, we also 
calculate the global metallicity [M/H] = [Fe/H] + 
Logj O.dSSfcy + 0.362) where fa = 10 (iSalaris et al. 
119931). Observed values of [a/Fe] from iGarretta et ah 
(1201 nal) are listed in Table |4] where available, as well as 
the values we assum e for our isochrone comparison (dis¬ 
cussed below in Sect. oil . We also list the corresponding 
global metallicity \M/H], inc l uding its uncertainty cal¬ 
culated following iNataf et al.l (j2013L see their eq. 7) and 
conservatively assuming a[a/Fe]=0A. In Fig. [6l we dis¬ 
play our fiducial sequences, shifted to the dereddened 
plane and color coded by [Fe/H]. 

4. MODELS VERSUS DATA ON THE RED GIANT BRANGH 
4.1. Adopted Models 

Our observations are compared to the following sets 
of evolutionary models, all of which have been updated 
recently, including the capability to generate ^-enhanced 
isochrones: 

1. Dartmouth St ellar Evolutionary Database (DSED; 
iDotter et al.ll2QQ8l ) isochrones with \a/Fe]=0.4:. 

2. Victoria-Regina (VR; IVandenBerg et ah] I2Q14D 
isochrones with [a/Fe]=0.4:. 

3. Bag of Ste llar Trac ks and Isochrones (BaSTI; 
IP let r infer ni et all I2QQ6I) canonical ^-enhanced 
isochrones. 

4. Yale-Yonsei (Y Y; lYi. Kim fc Demarauel I2QQ3I: 
iDemaraue et al.|[2QQ4D isochrones with [a/Fe] =0.4. 

5. Princeton-Goddard-PUC (PGPUG; 
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Fig. 5. — CMDs of our target clusters, with fiducial sequences overplotted in red and known variables as blue diamonds. In cases where 
a radial cut was necessary to reduce field contamination, the location of the radial cut in arcmin from the cluster center is given, and stars 
outside this radius are plotted in grey rather than black. Median color and magnitude errors in bins of 1 mag are shown along the right 
side of each CMD. 


IValcarce. Catelan fc SweigartI |2Q12[) isochrones, with 
[(a/Fe]=0.3 (this is the highest degree of a enhancement 
available over GGC parameter space). 

Any critical attempt to evaluate these models versus 
each other on equal footing is not truly an apples-to- 
apples comparison, as they all differ in various input as¬ 
sumptions (e.g. solar abundances, ^-enhanced chemical 
compositions, treatment of convection and atomic diffu¬ 
sion, color-Tg// relations). However, a direct comparison 


between different near-IR isochrones has yet to be pub¬ 
lished, and comparisons between individual isochrones 
and high-quality data are scarce. Here we take ad¬ 
vantage of our homogenous photometry of 12 GGCs 
covering a range of distance, reddening and metallic- 
ity to objectively explore how well existing isochrones 
fit observational data in the near-IR. That being said, 
in a few cases some modifications were made to the 
available default settings of various models in an at- 
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TABLE 4 

Cluster Parameters 


Gluster 

RA (J2000)'‘ 

Dec (J2000)'‘ 

(m — M)o^ 

E(B - V)'^ 

[Fe/H]‘= 

[a/Fe\^ 

observed 

[a/Fe] 

assumed 

[M/HY 

assumed 

A(V - lY 

NGG104 

00:24:05.71 

-72:04:52.7 

13.26 

0.023 

-0.76+0.02 

0.42 

0.4+0.1 

-0.47+0.08 

0.153+0.003 

NGC288 

00:52:45.24 

-26:34:57.4 

14.83 

0.013 

-1.32+0.02 

0.42 

0.4+0.1 

-1.03+0.08 

1.022+0.025 

NGG362 

01:03:14.26 

-70:50:55.6 

14.76 

0.023 

-1.30+0.04 

0.30 

0.4+0.1 

-1.01+0.09 

0.195+0.003 

NGG1261 

03:12:16.21 

-55:12:58.4 

16.08 

0.013 

-1.27+0.08 


0.4+0.1 

-0.98+0.11 

0.201+0.005 

NGG1851 

05:14:06.76 

-40:02:47.6 

15.42 

0.020 

-1.18+0.08 

0.38 

0.4+0.1 

-0.89+0.11 

0.234+0.011 

NGG2808 

09:12:03.10 

-64:51:48.6 

15.05 

0.183 

-1.18+0.04 

0.33 

0.4+0.1 

-0.89+0.09 

0.966+0.025 

NGG4833 

12:59:33.92 

-70:52:35.4 

14.19 

0.359 

-1.89+0.05 


0.4+0.1 

-1.60+0.10 

0.900+0.029 

NGG5927 

15:28:00.69 

-50:40:22.9 

14.57 

0.399 

-0.29+0.07 


0.2+0.1 

-0.15+0.10 

0.112+0.007 

NGG6304 

17:14:32.25 

-29:27:43.3 

13.99 

0.481 

-0.37+0.07 


0.2+0.1 

-0.23+0.10 

0.105+0.004 

NGG6496 

17:59:03.68 

-44:15:57.4 

14.95 

0.216 

-0.46+0.07 


0.2+0.1 

-0.32+0.10 

0.107+0.008 

NGG6584 

18:18:37.60 

-52:12:56.8 

15.71 

0.079 

-1.50+0.09 


0.4+0.1 

-1.21+0.12 

0.408+0.062 

NGG7099 

21:40:22.12 

-23:10:47.5 

14.72 

0.054 

-2.33+0.02 

0.37 

0.4+0.1 

-2.04+0.08 

0.872+0.006 


^ From IGoldsburv et al.l (|2010l) 

^ FromJD^ 

From^ _ 

^ From fCarretta et al.l (I2010al) 

® Calculated using [M/H] = [Fe/H] + Lo^f (0.638/a +0.362), where [Q;/Fe] = [a/Fe](assumed) 



Fig. 6. — Our fiducial sequences in the Mk, {J—Ks)o plane, color 
coded by [Fe/H] value. Median photometric errors in magnitude 
bins are shown along the right-hand side. 

tempt to maximize consistency. First, the near-IR col¬ 
ors and magnitudes output by the BaSTI and YY mod¬ 
els a re on the Johnson-Co usins-Glass photometric sys¬ 
tem (|Bessell fc BrettI [l988[) . so they were converted to 
the 2MASS pho tometric system using the relations of 
iCarpent^ (|2QQ1[ ). Second, the VR and PGPUC models, 
rather than assuming a AYjAZ relation (or having one 
optionally available), give the user the ability to interpo¬ 
late simultaneously in Y, [a/Fe] and [Fe/H] (or Z). In 
these cases, w e hav e chosen values of Y in order to be 
consistent with lDlOl increasing Y slightly with [Fe/H] to 
approximate a AY/AZ=1.4 relation. In any case, minor 
model-to-model differences in AY are inconsequential to 
our results since the location of the RGB in the near-IR 
is quite insensitive to He variations (see Sect. lOl) 

4.2. Cluster by Cluster Comparison 

In Fig. [71 we compare the five sets of isochrones to 
observed ISPI photometry for eight of our target clus¬ 
ters with the highest quality photometry spanning the 
observed metallicity range. For increased clarity, we 


show two panels for each cluster, where the left panel 
directly overplots the isochrones on the observed pho¬ 
tometry, and in the right panel we plot the color dif¬ 
ference with respect to the observed fiducial sequence 
as a function of M^. For this purpose, we conserva¬ 
tively approximate the uncertainty on the color of the 
fiducial sequence as the standard deviation of the color 
difference between the fiducial sequence and the ob¬ 
served RGB stars, calculated in moving bins of width 
0.5 mag in 7^5 using exclusively stars redward of the 
fiducial sequence to avoid contamination from the HB 
and AG Em. The resulting color uncertainty as a func¬ 
tion of magnitude is illustrated using dotted lines in the 
right-hand panel of Fig. [7] for each cluster. We employ 
models with [a/Fe]=0A (except in the case of PGPUG as 
noted above), but for the three most metal-rich clusters 
(NGG5927, NGG6304, NGG6496, -0A6<[Fe/H] <0.29) 
we use models wit h [a/ Fe]=0.2. This choice is supported 
by both ETaandm and is consistent with observed 
trends of [a/Fe] vs. [Fe/H] a s a function of Gala cto- 
centric radius and scale height (jHavden et al.l l2015l and 
references therein) and the ’’disklik e” enrichment sce¬ 
nario employed in earlier studies (e.g. lFerraro et al.l! 19991: 
[Ferraro. Valenti fc Origli^l2006F 

Total agreement, either between models and data or 
among various sets of models, would be surprising due in 
part to the v arious model-to-model differences mentioned 
in Sect. EH However, it appears that the DSED and VR 
models reasonably reproduce the observed morphology 
of the observed RGBs, although deviations in morphol¬ 
ogy from the fiducial sequences are still seen (as with 
all models investigated here). In addition, these devia¬ 
tions appear to become more drastic at low {[Fe/H] <-2) 
metalliciti es. While this result is fully consistent with the 
findings of iBrassenr et al.l (|2Q1Q[) using optical-IR colors, 
we caution that at least in our case, the formal statis¬ 
tical significance of these morphological deviations may 
not be high when photometric errors as well as calibra¬ 
tion uncertainties are taken into account. 

In cases where few stars are available (typically close to the 
RGB tip), the size of the magnitude bin was expanded until at 
least 5 stars were available to calculate the standard deviation 
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Fig. 7. — Comparison between isochrones from five sets of evolutionary models and observed cluster photometry. For each cluster, the left 
panel illustrates the isochrones overplotted on the photometry, with median photometric errors given along the right hand side. The right 
panel illustrates the color difference between the models and the observed fiducial sequence as a function of absolute Ks magnitude, where 
the vertical solid line represents equality and the dotted lines indicate the estimated color uncertainty of the fiducial sequence. Models are 
color coded as indicated in the right hand panel for each cluster. 

To quantify the color difference A(J — Ks) between the observed RGB fiducial sequence and the five mod- 
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els shown in Fig. [7| for each cluster, we have calcu¬ 
lated the mean and standard deviation (weighted using 
the observed errors shown as dotted lines in Fig. [7]) of 
A{J — Ks) in evenly spaced magnitude bins over the 
cluster RGBs (Mx<0). In this way, the mean quantifies 
the extent to which a model matches observed colors^ 
averaged over the RGB, while standard deviation gauges 
the match in terms of GMD morphology. These statis¬ 
tics are given for all target clusters in Table [5l and are 
summarized across all clusters using the median and the 
median absolute deviation (MAD) in the last two rows. 
The results indicate that the DSED and VR models do 
the best job of reproducing the observed GMD morphol¬ 
ogy, although BaSTI and PGPUC do nearly as well in 
this sense, and the BaSTI models actually best agree 
with the observed fiducial colors in the mean. 

We di scus s practical ramifications of this result further 
in Sect. 16.2| and now explore the effects of variations on 
[a/Fe] and helium predicted by the models, in order to 
gain insight into the role of uncertainties in these quan¬ 
tities. 

4.3. Model Predietions: [a/Fe] Variations 

The difference between scaled solar isochrones versus 
those which are ^-enhanced (at fixed Y and [Fe/H]) are 
illustrated in Fig. [HI There, the {J—Ks)o color difference 
between isochrones with [a/Fe]=0.4 and 0 are plotted in 
the sense (^-enhanced—scaled solar), for a regular grid of 
[Fe/H] values spanning the values of the target clusters. 
To illustrate the influence of a-enhancement in the near- 
IR on the main sequence, we extend these plots faintward 
beyond the lower main sequence knee (MSK). A couple of 
features are evident: First, all of the models except for 
BaSTI give schematically similar predictions regarding 
the effect of a-enhancement on the upper RGB at low to 
intermediate {-2<[Fe/H] <-l) metallicities. Specifically, 
an enhancement of [a/Fe]=0.4 dex shifts the upper RGB 
redward by 0.05 mag or more in {J — Ks) color (depend¬ 
ing on the model and the [Fe/H] value), larger than seen 
anywhere else in the GMD with the exception of the lower 
{Mk>^) main sequence. Meanwhile, the lower RGB 
(0<Mx<2) remains less affected save for a slight ^0.02 
mag redward shift. However, approaching solar metallic- 
ity, the model predictions diverge somewhat, with some 
(DSED, VR) predicting a negligible effect on the RGB 
tip but maintaining a redward shift over the remainder 
of the RGB. 

4.4. Model Predietions: Helium Enhaneement 

Much effort has recently been devoted to the issue of 
a spread or enhancement in helium among GGG (sub- 
) populations, as there are important implications for 
GGG formation. Therefore we briefly examine to what 
extent a change in helium abundance influences near-IR 
GGG CMDs. The result of an increase of AV=0.04 (at 
fixed [a/Fe] =0.4) is shown in Eig. |9l again at a range 
of [Fe/H] values. Only the VR models are shown, al¬ 
though the results are nearly identical for other models, 
which is that near-IR colors on the RGB are minimally 
affected (A(J —i^ 5 ')< 0.01 at fixed M^), with essentially 
no consequences for the RGB morphology. As it has been 
demonstrated that the MST O and MSK can b e used as 
age indicators in the near-IR (jBono et al.lDOlOl ) , we point 



Fig. 8.— The color difference between an a-enhanced isochrone 
and a scaled solar isochrone for each model as a function of ab¬ 
solute Ks magnitude. The difference is plotted in the sense {a- 
enhanced—scaled solar), and color coded by model as in previous 
figures. 

out that while helium enhancement shifts the MSTO and 
the RGB very slightly blueward in color, these shifts are 
equivalent to <0.01 mag. Similarly, the magnitude of the 
MSTO and MSK are shifted faintward by nearly equal 
amounts (<0.1 mag, see Eig. [10]). Eor this reason, a 
modest helium spread or enhancement in the near-IR is 
of little consequence for age d eterminations employing 
relative GMD indices fe.g. IVlSl l. 


4.5. Photometrie Indiees 


Previous studies ( 

Erogel. Gohen & Perssod [19831 

Gohen & SleeDer[[1995[: 

Eerraro et al.[ 20001 [Valenti et al.l 

2004aD. cognizant of t 

he fact that the RGB becomes 
to metallicity at higher lumi- 
set of photometric indices to 
RGB as a function of metallicity 

increasingly sensitive 
nosities, developed a 
characterize the upper 


in the near-IR. These indices include (dereddened) color 
at fixed (absolute) magnitude, specifically (J — Ks)o 
at Mx=(-3,-4,-5,-5.5), as well as absolute magni - 
tude Mk at (J — i^s)o=0.7. As [Eerraro et al.[ (|2QQ0[ ) 
point out, the use of indices constructed purely using 
IR filters has the advantage of decreased sensitiv¬ 
ity to uncertainties in both distance, due to the 
near-v ertical slope of the RGB (e.g. [Valenti et al.[ 
2004aD. and reddening (E(J — Ks)=0.^33E(B — V); 
Gasagrande fc VandenBei^ 120141) . These indices can 
be employed to meas ure cluster distances, reddening s 
and metallicities (e.g. [Eerraro. Valenti fc Origlig][2006[) . 
and also provide another approach to quantifying 
discrepancies between isochrones and observed fiducial 
sequences. Therefore, in Eig.[TT]we present new relations 
between these photometric indices as function of both 
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TABLE 5 

Color Difference A(J - Ks) Between Fiducial Sequences and Models 


Cluster 

DSED 

Mean a 

VR 

Mean a 

BaSTI 

Mean a 

YY 

Mean a 

PGPUC 

Mean a 

NGC5927 

-0.049 

0.018 

-0.071 

0.009 

-0.036 

0.016 

0.014 

0.033 

-0.053 

0.028 

NGC6304 

-0.021 

0.017 

-0.040 

0.017 

-0.014 

0.022 

0.034 

0.027 

-0.044 

0.014 

NGC6496 

0.002 

0.016 

-0.017 

0.019 

-0.004 

0.026 

0.039 

0.044 

-0.023 

0.023 

NGC0104 

-0.001 

0.011 

-0.012 

0.011 

0.023 

0.016 

0.046 

0.028 

-0.012 

0.012 

NGC1851 

-0.023 

0.008 

-0.031 

0.007 

-0.003 

0.020 

0.031 

0.040 

-0.035 

0.006 

NGC2808 

-0.008 

0.019 

-0.014 

0.008 

0.020 

0.012 

0.059 

0.033 

-0.026 

0.021 

NGC1261 

-0.017 

0.017 

-0.032 

0.016 

-0.007 

0.019 

0.026 

0.035 

-0.033 

0.020 

NGC0362 

-0.005 

0.020 

-0.020 

0.009 

0.004 

0.010 

0.037 

0.026 

-0.021 

0.020 

NGC0288 

-0.012 

0.008 

-0.033 

0.021 

-0.008 

0.033 

0.016 

0.054 

-0.032 

0.012 

NGC6584 

-0.008 

0.018 

-0.028 

0.013 

-0.000 

0.018 

0.038 

0.036 

-0.029 

0.019 

NGC4833 

0.006 

0.007 

-0.016 

0.008 

0.013 

0.013 

0.057 

0.037 

-0.015 

0.009 

NGC7099 

-0.016 

0.011 

-0.046 

0.012 

-0.022 

0.014 

0.008 

0.030 



Median 

-0.008 

0.017 

-0.028 

0.012 

-0.003 

0.018 

0.037 

0.035 

-0.029 

0.019 

MAD 

0.009 

0.004 

0.012 

0.003 

0.011 

0.005 

0.011 

0.005 

0.006 

0.005 



-0.05 0.00 0.05 -0.05 0.00 0.05 -0.05 0.00 0.05 

A(J-Ks) 


Fig. 9.— The J — Ks color difference predicted for a helium 
enhancement of Ay=0.04 from [a/Fe\=^A VR models at a range 
of [Fe/H]. The difference is plotted in the sense (helium enhanced- 
helium normal). 

cluster [Fe/H] and global metallicity [M/H] as well as 
a comparison with the five sets of evolutionary models 
listed in Sect. 14.11 The uncertainties i n de reddened 
color are calculated as described in Sect. 14.21 and our 
linear fits are performed taking uncertainties in both 
axes into accounlF^. The resulting equations and the 
rms deviation of the residuals is given in each panel of 
Fig. m as well as inverted versions of these equations 
(with [Fe/H] or [M/H] as the dependent variable) 
and the corr e spondi ng x-axis rms. The relations of 
1 Valenti et al.l (l2nn4aL transformed to the metallicity 
scale of lCQ9l . are overplotted as dotted lines for com- 

see http://idlastro.gsfc.nasa.gov/ftp/pro/math/fitexy.pro 



Fig. 10.— The effect of helium enhancement in the near-IR, using 
VR isochrones with a fixed [Q;/Fe]=0.4 and [Fe /shown 
for y=0.25 (solid line) and y=0.29 (dotted line). 


parison. It appears that a linear fit is a reasonable 
representation of {J — Ks)o color as a function of [Fe/H] 
and [M/H], with a slope that increases wi th luminosity, 
in general accord with p revious results (jFerraro et al l 
I2QQQI: IValenti et al.ll2QQ4al) . 

Dereddened color at fixed absolute magnitude can also 
be used to characterize the RGB, as these indices sam¬ 
ple a luminosity range which is dependent on metallicity 
(this is apparent, for example, in Fig. [6]). In Fig. [12] we 
present linear fits of Mk at fixed (J — i^5')o=0.7, again 
as a function of [Fe/H] and the global metallcity [M/H]. 
In this case, the observational magnitude uncertainty has 
been calculated by multiplying the observed color uncer¬ 
tainty by the slope of the fiducial sequence at the cor¬ 
responding magnitude. The resulting fits are shown in 
Fig. [T^ again with [ a/Fe] =^A and 0.0 models and the 
fits of lValenti et al.l (|2Q04a[ ) shown for comparison. 

Indices of color at fixed magnitude, unlike magnitude 
at fixed color, are extremely sensitive to cluster redden¬ 
ing values, as even the most metal- rich GGCs have RGB 
slopes of \5{J — Ks)/SKs\<0.12 (| Valenti et al.l l^04ai 
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Fig. 11.— Indices of dereddened (J — Ks)o color at fixed absolute magnitude Mk=- 5.5, -5, -4 and -3 (top to bottom) as a function of 
[Fe/H] (left) and [M/H] (right). Vertical error bars represent measurement uncertainties only and do not include uncertainties in cluster 
distance and reddening (see text for details). Predictions of evolutionary models are overplotted using the same colors as in Fig. [3 where 
models with [Q;/F'e]=0.4 are shown as solid lines and models with [a/Fe]=0 are shown as dashed lines. A linear fit is s hown a s a black solid 
line , wit h the coefficients and the rms deviation given in the upper left of each panel. The relations of I Valenti et al.l (I2004alb converted to 
the IC09I metallicity scale, are overplott ed a s grey lines in the left panels. The open diamond represents the value for NGC 5927 assuming 
a metallicity of [Fe/H]=-0.5 (see Sect. 16.2|) . 


I2Q1QI: [Chun et al.1l2QlQ[) . This is an advantage when us¬ 
ing color at fixed magnitude to describe the RGB, since a 
moderate uncertainty in [Fe/H] or [M/H] does not have 
drastic consequences for the determination of distance 
and reddening. For example, the slope of the relation in 
the upper left panel of Fig. [TT] implies that a[M/H]=0.3 
dex translates to aAx ^0.03. On the other hand, an 
uncertainty of aE{J — Ks) ^0.02, or aE{B — V) ^0.04, 
could account for a scatter of ^0.2 mag in the relations 
presented in Fig. [121 although in the present case obser¬ 
vational uncertainties appear to dominate. Interestingly, 
the quality of the linear fits in Figs. [TT] and [12] improve 
at higher luminosities when the global metallicity [M/H] 
rather than [Fe/H] is employed, demonstrating that the 
upper RGB in the near-IR is quantifiably affected by 
cluster [a/Fe] consistent with the model predictions in 


Sect. IJL Our linear fit s generally compare well to those 
of I Valenti et al.l (j2QQ4a[ ) given their uncertainties, and we 
find nearly identical slopes in the case of color at fixed 
magnitude (Fig. [TT]) . While zero point offsets are seen, 
the se offsets are not lar ge compared to the residuals of 
the I Valenti et al.l (j2QQ4al) fits, and can likely be explained 
by current uncertainties in cl uste r distances and redden¬ 
ings (discussed later in Sect. 16.1]) . 

5. HORIZONTAL BRANCH AND RED GIANT BRANCH 
BUMP MAGNITUDES 

5.1. Near-IR Horizontal Branch Magnitude 

In order to measure the magnitudes of our target clus¬ 
ter HBs, we first divide our sample into ’’red HB” (RHB) 

and ’’blue HB” (BHB) clusters^ _ This division is made 

using the A(V — I) index of IDIQI (given in Table [4]), 
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Absolute magnitude Mk at fixed color (J — Ks)q=^-7 as a function of [Fe/H] (left) and [M/H] (right). Symbols are as in 


Fig. 12.— 

Fig. El 

which is less vulnerable to saturat ion at its extremes than 
{B-R)/{B + V+R) (e.g.fig. 2of|Dia). We consider the 
seven clusters with (5(1/ — /)<0.3 as RHB clusters, and 
the remainder as BHB clusters, for the straightforward 
reason that all of the RHB clusters have HB LFs which 
show reliably detected peaks in magnitude whereas the 
BHB clusters have HBs which do not create a significant 
peak in their luminosity functions. The HB morphol¬ 
ogy of the BHB clusters is therefore not easily charac¬ 
terized using magnitude alone due to the diagonality of 
BHBs in near-IR CMDs (for this reason, we have in¬ 
cluded NGC 6584 among the BHB clusters despite its 
value of 5{V — /)=0.408) so we present HB magnitudes 
for only the RHB clusters. 

The location of the HB in J and Ks magnitude is quan¬ 
tified by constructing a luminosity function (LF) from 
the observed cluster CMD, and then identifying the loca¬ 
tion of the peak in this LF. To mitigate the effects of bin¬ 
ning on the final LF, we construct a multi-bin histogram 
that is the average of 10 individual LFs in which the bin 
starting points were shifted by 0.1 times the chosen bin- 
size and these 10 LFs are averaged. For clusters which 
have at least a portion of their comparison fields out¬ 
side of their Harris (1996, 2011 revision) tidal radii (in¬ 
cluding the most heavily contaminated GGCs NGG5927, 
NGG6304 and NGG6496) this procedure was repeated 
for the comparison field CMD, before scaling the field 
LF by the relative cluster to field area and subtracting it 
from the cluster LF. An example is illustrated in Fig. [131 
The uncertainty in the location of the LF peak is quan¬ 
tified using a simple Monte Carlo procedure wherein for 
each of 1000 iterations, the magnitude of each star in 
each filter is offset by a random amount drawn from a 
Gaussian distribution with a standard deviation equal 
to the photometric error. The construction of the LFs is 
then repeated, including the generation of the multi-bin 
histograms, and the location of the peak magnitude in 
each iteration reported. The uncertainty in the magni¬ 
tude of the LF peak is then the standard deviation of the 
reported LF peak locations over the 1000 Monte Carlo 
iterations. Typically, this uncertainty was smaller than 
a resolution element of the LF (0.01-0.03 mag), so to be 
conservative, the two values were added in quadrature 



Mk 

Fig. 13. — An example luminosity function for the case of NGC 
6304, the most heavily contaminated cluster in our sample, shown 
on a linear (top) and logarithmic (bottom) scale. After construct¬ 
ing a multi-bin LF from the cluster CMD (dark grey line), an LF is 
constructed for a comparison field, scaled by the ratio of cluster to 
field area (light grey line) and subtracted to obtain the decontam¬ 
inated cluster LF (black line). The observed LF peak is indicated 
by a vertical arrow in both panels. 

to obtain the final HB magnitude uncertainty, given in 
Table |6l 

We plot the reported HB magnitudes as a function 
of [Mjm in Fig. (Ill color coded by the cluster age 
from Id id In acidition, we overplot the predictions of 
iSalaris fc Girardil (j2QQ2f ) as dashed lines aft er convert¬ 
ing t heir Mk values to the 2MASS system (jCarpenterl 
Ml), since this remains the only study directly predict¬ 
ing HB magnitudes in the near-IR as a function of cluster 
age and metallicity which extends to the (low metallic- 
ity, high age) parameter space occupied by GGGs. At 
the metallicity of 47 Tuc and higher, ou r results are 
in agreement with ISalaris fc Girardil (j2QQ2l ) given uncer¬ 
tainties on a bsolu te age (>0.5 Gyr from measurement 
errors alone; Ima), but at lower metallicities, the mod- 
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Fig. 14. — Horizontal branch absolute magnitude in the Ks (top) 
and J (bottom) filters versus cluster [M/if] for RHB cluster s. In 
the upper panel, the predictions of ISalaris fc Girardil (|2002li are 
overplotted as dash ed lines after tra nsformation to the 2MASS 
photometric system (|CarDenterll200lT) . and cluster and model pre¬ 
dictions are both color coded by age as indicated in the top left of 
the upper panel. 

els overestimate the iC-band luminosity of the HB by as 
much as 0.1 mag in the case of NGC 1851. This issue 
should clearly be revisited with a larger sample of GGCs, 
and inter- an d intra-cluster abundanc e variations not in¬ 
cluded in the ISalaris fc Girardil (j2QQ2[ ) models could pos¬ 
sibly account for this discrepancy, but perhaps a more 
likely possibility is simply the remaining uncertainty in 
the GGC distance and age scales (this is discussed fur¬ 
ther in Sect. lOTi . However, a dependence of Mk{HB) 
on metallici ty for GGGs has been observed at least since 
the study of iGrocholski fc Saraiedinil (j2QQ2[ ) , and the ob¬ 
served HB magnit udes are not the culprit: The m edian 
HB magnitudes of iGrocholski fc Saraiedinil (j2QQ2[ ). mea¬ 
sured directly from 2MASS, differ from ours by <0.04 
mag (after accounting for t he offset of 0.044 mag be¬ 
tween photometric systems; iCaroent^ l200lD for both 
GGGs included in their study (47 Tuc and NGG 362). 

We _ note that ou r c omparison to 

iGrocholski &: Saraiedinil (j2002[ ) and ISalaris fc Girardil 
(j2002f ) is not strictly homogenous, as the former use the 
median of a CMD-selected region to characterize the HB 
magnitude and the latter give a mean HB magnitude. 
However, we found that employing any of these alternate 
techniques yielded HB magnitudes which differed from 
the LF peak by ^0.02 mag in the mean, so that in 
this case, the choice of methodology for quantifying the 
HB peak does not significantly impact our conclusions, 
although radial population gradients which vary with 
HB color (fo r example due to lig ht element abundance 
variations cf. iGratton et al.ll^oTol ) could play a role. 

5.2. The Red Giant Braneh Bump (RGBB) 

The magnitude of the RGBB is m e asure d using a pro¬ 
cedure similar to that of iNataf et al.l (j2Q13f ) . Specifically, 
the RGB LF of a cluster is built using stars selected from 
the CMD (e.g. see Fig. [15]) . and multibin histograms are 



0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 1.2 

(J-Ks) (J-Ks) 

Fig. 15.— Examples of the selection of RGB stars for the con¬ 
struction of the RGB LF, shown for the cases of NGG 104 (47 Tuc; 
left) and NGG 7099 (M30; right). All cluster sources are shown in 
grey, and stars used to construct RGB LFs are shown in black for 
each cluster. 

used to construct the LF as described in the previous sec¬ 
tion. Again, the field LF (constructed from an identical 
GMD region as the cluster LF) is subtracted where fea¬ 
sible. Next, an exponential plus Gaus sian function is fit 
to the resulting cluster RGB LF (e.g. iNataf et ani2Q13D 
to measure the location of the RGBB, and an example of 
this procedure is shown in Fig. [161 Similar to Sect. 15.11 
Monte Carlo simulations (over 1000 iterations) are used 
to quantify the uncertainty on the RGBB magnitude in 
each bandpass, wherein for each iteration, each star is 
offset by a Gaussian deviate of its photometric error, the 
multibin LF is reconstructed, an exponential plus Gaus¬ 
sian function is refit, and the resulting best fit location of 
the Gaussian peak is reported. The errors reported here 
were then calculated as the quadrature sum of the fit 
uncertainty to the observed RGBB plus the standard de¬ 
viation of the best fit RGBB magnitude across the Monte 
Carlo iterations. 

The resulting bump magnitudes and their uncertain¬ 
ties are listed in Tabled and the RGBB absolute magni¬ 
tude in each filter, Mj{RGBB) and Mk{RGBB) are 
shown in Fig. [T7| a s a function of both [Fe/H] and 
[M/H]. Similarly to lValenti et al.l (|2QQ4b[) . we were un¬ 
able to detect the RGBB in NGC 7099. However, by 
matching our ne a r-IR d ata to the optical photometry by 
iSaraiedini et al.l (|2QQ7[) we were able to combine t he V- 
band magnitude of the RGBB from INataf et al.l (|2Q13[ ) 
with its (V — Ks) color to calculate the Ks magni¬ 
tude of the RGBB, shown in Fig. EL Next, w e com ¬ 
bine d our sample with cl usters from iCho fc Led (12001) 
and I Valenti et al.l (|2004b[ ) which have high-qualit y dis¬ 
tances, reddenings and spe ctroscopic metallicities (|C09l: 
iCarretta et al.lHoiOal: IDIOF and performed a quadratic 
fit for direct compar i son with p revious results. The rela¬ 
tion of IValenti et al.l (|2004b[fn converted to the lC09l scale 
appears to deviate from ours at high (near solar) metal- 

The coefficients given by I Valenti et al.l (I2004bl) for the 
Mk(RGBB) — [Fe/H] relation differ slightly between their fig. 2 
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Fig. 16.— Example iF^-band RGB LF for NGC 6304, shown on 
a li near (top) and logarithmic (bottom) scale. Symbols are as in 
Fig. I13l and the exponential plus Gaussian fit used to quantify the 
RGBB magnitude (indicated by the vertical arrow) is shown as a 
dotted line. 

licity, and we urge further multi-object spectroscopic 
studies of the most metal-rich G GCs to confirm this re - 
sult. However, transforming the IValeuti et al.l (l2QQ4bD 
relatio n on the met allicit y scale of iGarretta fc GrattonI 
(|1997f ) to the scale of lCQ9l comes with the caveat that the 
transformation between these scales is in fact poorly de¬ 
fined at high metallicities, as none of the 13 clusters em¬ 
ployed to derive the linear transformation bet ween scales 
are more metal rich than 47 Tuc (see fig. 8 ofEHl). 

6. DISGUSSION 

6.1. Uncertainties in Cluster Parameters 

To maximize homogeneity we consistently use the lat¬ 
est compilation of GGC distances and reddenings from 
iDlQl . which are based on fits of DSED isochrones to deep 
optical photometry from iSarajedini et al.l (j2QQ7f ). While 
a detailed discussion of the GGC distance scale is be¬ 
yond the scope of this investigation, recent studies esti- 
mate current uncertainti es at the level of ^0.10-0.15 mag 
(|V13l : IVandenBer jl2Q13L and references therein), consis¬ 
tent with a dir ect comparison between the distances re¬ 
ported bv Id id versus those o btained via subdwarf fittin g 
to the GGC main sequences (jCohen fc Saraiedinill2Ql^ . 

In the current co ntext, we revisit the comparison from 
iCohen et al.l (l2Ql4 see their fig. 12) between the results 
of ima and IV13l with regard to GGC distances, red¬ 
denings, and metallicities. This comparison is especially 
useful since both inves tigations employed iden tical pho¬ 
tometric catalogs from ISarajedini et al.l (j2QQ7[ ) , but dif¬ 
ferent str ategie s to obtain cluster parameters. On the 
one hand, IDIQI optimized their isochrone fits for the un¬ 
evolved main sequence and the lower RGB, using values 
from the iHarrisI ()1996f ) catalog as initial guesses and al¬ 
lowing slight variations to [Fe/H]^ {m — M)Fd>iAW and 
E{F^ t}UW — F8141T) when necessary. On the other 
hand, IV13l fit model ZAHBs to the observed lower en- 

and their eq. 3, and the former are used in Fig. 1171 


velope of clus ter H Bs to measure ages, restricting these 
models to the lCQ9l [Fe/H] values and a two part linear 
[a/Fe]-[Fe/H] relation. 

In Fig. [iHl we present a comparison between the 
(m — M)o, E {B — V), [Fe/H] and age values employed 
by Ima and MB- The distance discrepancy between 
the two studies is generally within the aforementioned 
uncertainties (^0.1 mag), and the reddening values are 
in particularly excellent agreement. We also investigate 
the relevance of the difference in m etalli cities between 
the spectroscopic values report ed by lCQ9l which we em¬ 
ploy and the values used by iDlOi for their isochrone 
fitting. The two sets of valu es ar e generally in good 
agreement (see Fig. HHI, and Ima discuss the role of 
uncertainties in [Fe/H] in some detail. Ho wever, in 
two cases the metallicities employed by Ima for their 
isochrone fits differ signi ficantly from fCQ9l : F irst, for 
NGC 5927 ([FeJH] Scm =-0.29+()m: [Fe/H]^^=- 
0.5), the most metal rich cluster in our sample, dis¬ 
tance and reddening are unlikely to play a role as both 
E{J — Ks ) an d (m — M)k agree to <0.01 mag be¬ 
tween Ima and IV13l . Furthermore, transformations of 
earli er [Fe/H] measurements for this cluster to the 
C09l scale rest fairly heavily on the value obtained by 
Garrett a et al (j2001l ) for NGG 6528 (see figs. A.l and 
A.2 of iCOOl ). while subsequent observations of NGG 6528 
imply lower val ues and a substantial (^0.2 dex) spread 
in [Fe/H] fe.g. IZoccali et~al] l2004 iSobeck et al.l 120061 : 
iManro et al.l l20 1^ . The ot her si gnifi cant discrepancy in 
cluster metallicity b etwee n lC09l and iDlOl is the case of 
NGG 4833 ([Fe/i7](|C09|)=-1.89±0.05; [Te/i7] (lmol) =- 
2.4). Recent spectroscopic evidence cl arifies the sit¬ 
uation to some degree: On one hand, IGarretta et al.l 
(j2014l ) repo rt [F e/i7]=-2. Q15T0.09 dex from UVES spec- 
tra on their lC09l scale, but lR.oederer fc ThompsonI (j2015f ) 
find [Fe/i7]=-2.25±0.02 from Fe I lines and [Fe/H]=- 
2.19±0.013 from Fe H lines, although they show that the 
use of different model atmosphere grids and line analysis 
codes can fully account for the discre pancy between their 
results and the IGarretta et all (j2014j ) value. Fortunately, 
the aforementioned issues are of little consequence for the 
results of Figs. [HI [12] and [T3 In fact, if we use values 
of (m — M)o a nd E {B — V) from IV13l and/or [Fe/H] 
values from the iDlQl isoc hrone fits (rather than the spec¬ 
troscopic values of ICQ9[ ). the fits in Figs. [TTl \V2\ and [TTI 
are unaffected to within their rms deviations. Alterna¬ 
tively, we could have simply excluded NGC 4833 and/or 
NGC 5927 from our fits, but this also turns out not to 
impact the fits beyond their uncertain ties. _ 

Our result from Sect. I01 that the iSalaris fc Girardil 
(j2QQ2[ ) models overestimate the HB luminosity for metal- 
intermediate {[M/H] <-0.8) RHB clusters is also ro bust 
to the choice of di stances, reddenings, and ages from I VI 3l 
rather than IDIQI . This is due at least parti ally to the 
fact that although the abs olute ages of IV13l are gener¬ 
ally younger than those of IDIQI . the relative ages of our 
target clusters are unchanged within their quoted un¬ 
certainties regardless of which set of absolute ages one 
assumes (see Fig. [TH] bottom panel). In f act, the use 
of di stances, reddenings and ages from rather than 
IDIQI actually increases the discrepancy between observed 
and predicted Mk at lower metallicities: Taking reported 
age and metallicity uncertainties into account, the shift 
required to bring observed Mk{HB) values into accord 
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Fig. 17.— Red giant branch bump absolute magnitude in the J (upper panels) and Ks (lower panels) filters versus cluster metallicity, in 
terms of \Fe/H] (left) and \M / H] (right). Filled black circles represent values for ISPI target clusters, grey circles are from I Valenti et 
(|2004bl b and black open diamonds are from ICho &: Led (|2002l b The black square represents the RGBB location of NGC7099 predicted from 
optical-infrared photometry (see text) and has not been included in the fits. The sol id line is a quadratic fit (with the resultin g coe fficients 
given in the upper right of each panel), and the dotted line represents the relation of I Valenti et al.l (I2004blb transformed to the TG09l [Fe/H] 
scale. 


TABLE 6 

Red Giant Branch Bump and Horizontal Branch Magnitudes 


Gluster 

M(RGBB) 

M{HB) 


J 

Ks 

J 

Ks 

NGG0104 

12.728T0.010 

12.072T0.009 

12.492T0.013 

11.979T0.013 

NGG0288 

13.736T0.035 

13.150T0.034 



NGG0362 

13.754T0.026 

13.136T0.029 

14.155T0.011 

13.680T0.011 

NGG1261 

14.997T0.009 

14.420T0.009 

15.418T0.022 

14.928T0.026 

NGG1851 

14.424T0.025 

13.808T0.028 

14.775T0.012 

14.340T0.011 

NGG2808 

14.253T0.013 

13.530T0.012 



NGG4833 

12.917T0.014 

12.133T0.028 



NGG5927 

14.585T0.027 

13.760T0.038 

14.036T0.013 

13.254T0.011 

NGG6304 

14.211T0.017 

13.295T0.029 

13.577T0.022 

12.709T0.024 

NGG6496 

14.783T0.024 

14.015T0.035 

14.288T0.036 

13.666T0.030 

NGG6584 

14.624T0.011 

13.999T0.010 




with the lSalaris fc Girardil (j2QQ2[ ) models at the Icr level 
is AMx(i^5)=-0.15 and -0.17 mag for NGC 362 and 
NGC 1851 respectively, compared to AM/u(i7^)= -Q.Q5 
and -0.10 mag respectively (see Fig. using the Ima 
distance, reddening and age. Howeve r, we recall that 
the distances and reddenings given by IV13l were deter¬ 
mined by assuming the lower envelope of the observed 
HB as the ZAHB (at the high-mass end). This proce¬ 
dure comes with the risk that any dependence on metal- 
licity of the magnitude difference between the model 
ZAHB and the observed lower envel ope of the HB (as 
disc ussed extensiv ely, for example, bv lFerraro et al .111 9991 
and iGatelanI 120091 ) could systematically bias ou r resu lts. 
For this reason, the distances and reddenings ofima ap¬ 
pear a somewhat more objective set of values with which 
to investigate a relation between HB absolute magni¬ 
tude and metallicity, in the sense that their isochrone 
fits, optimized for the cluster main sequences and lower 


RGB, give distances and reddenings perhaps more “in¬ 
dependently” of the HB0. We also caution that the in¬ 
terplay between HB morphology and other cluster pa¬ 
rameters, including age, metallicity, light element abun¬ 
dances, and helium abundance, appears e xtremely com¬ 
plex and is under ac t ive investigation fe.g. [Catelanl[2009l: 
iGratton et al.l I2QTqI : IDIQI : iMilone et al.l 1201^ . As our 
data are affected by central incompleteness, we cannot 
exclude the possi bility that radial gradients in HB (sub-) 
populations (e.g.lNataf et al.ll2QlirKrQgsrnd et al.ll2Q13E 

IVanderbeke et aill2Q15f ) play a role. Therefore, a more 

detailed comparison between observed and synthetic HBs 
is a topic better tackled using high spatial resolution 
imaging of cluster cores. 

Althoug h maybe not in the most s trict sense: The distance 
moduli in the lHarrii (I1996l i catalog which lPTol use as initial guesses 
are based on an empirical Mv{HB) — [Fe/H] relation. 
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[Fe/H] (C09) 


Fig. 18. — Comparison of difference in (m — M)o, E{B — V), 
\Fe/ H] an d cl uster age in Gyr (top to bottom resp ectively) be¬ 
tween |DT0 and lVl^ as a function of [Fe/H] from lCO^ . The dashed 
line in each plot represents equality, and our target clusters are 
overplotted in black. 


Differences in cluster distances and reddenings may 
also be respo nsible for minor differ ences between our fits 
and those of I Valenti et alJ (j2QQ4allbl ) in Figs. nmn] and 
[m In Fig. [191 we plot t he di fference in (m — M) n and 
E{B — V) between IDIQI and iFerraro et al.l (1199911. t he 
latter of which was employed by IValenti et al.l (|^Q4allbl ) 
to calibrate their near-IR relations. Again, our target 
clusters are shc )wn as filled black ci rcles, and we have 
Valenti et al.l 


overplotted the 
ters common to 


ison 


between ID 


, _ _ _ calibrating clus- 

DlQlwith d iamonds. Unlike the compar- 


3 and I VI 31 . variations in cluster distance 


and reddening range up to >0.1 mag in E{B — V) and 
>0.2 mag in (m —M)o (albeit in the opposite sense: clus¬ 
ters with longer distance s have smaller redden ing values). 
The lower reddenings of iFerraro et al.l (Il999[) are a likely 


contrib utor to the blueward offset of the IValenti et al.l 
(|2004af ) relations in Fig. [TT] (recall that near-IR color at 
fixed magnitude is more sensitive to reddening than dis¬ 
tance) , as well as the brightward offset of their relation in 
Fig. [T2l However, it is possible to explain the offset b e- 
tween our relations and those of IValenti et al.l (j2004allbl ) 
without invoking differences in calibrating cluster dis¬ 
tances and reddenings given the calibration uncertain¬ 
ties, measur ement uncertainties a nd fit quality obtained 
here and by IValenti et al.l 


6.2. Implications for the Use of Isochrones 

In Sect. im we found that while models are generally 
offset redward from our fiducial sequences in {J — Ks) 
color, the VR and DSED models were most success¬ 
ful at reproducing the morphology of the RGB. In fact. 



Fig. 19.— Comparis on of differenc e in (rn — M )p (to p) and F(B — 
V) (bottom) b etwee n IFerraro et al.l (|1999lf a ndlD10| as a function 
of [ Fe/H ] from lCO^ . Clusters present in both IFerraro et HI (|1999lf 
and IDIOI are shown as grey filled circles, our target clusters are 
qverplotted as filled black circles, and calibrating clusters from 
IValenti et al.l (|2004ai rbl) with distances and reddenings given bv lV13l 
are overplotted with black diamonds. The horizontal dotted line 
represents equality. Note the difference in y-axis scales as compared 
to Fig. [18] 


if we simply apply a fixed color offset to the DSED 
and VR isochrones, they reproduce the observed clus¬ 
ter sequences down to the main sequence for metal- 
intermediate clusters, as shown in Eig. [20] 

There, we have applied the mean color offset from Ta¬ 
ble [5] to the DSED and VR isochrones to correct for the 
difference between the observed color and that given by 
the isochrone. This mean color offset and its standard 
deviation (calculated as described in Sect. 14.2p is given 
for each cluster in each panel of Eig. [20l and is plot¬ 
ted as a function of cluster [Ee/H] in Eig. [2TJ In the 
case of the VR models, the blueward shifts which we 
find necessary were also found at optical wavelengths by 
IV13l in the sense that their ZAHB fits yielded MSTO 
and RGB colors which were too red, but these shifts, 
which were all <0.03 mag in {E606W — F814IU), are 
too small alone to explain the present color discrep¬ 
ancies of A{J — Ks){VR) ^0.04, which translates to 
A(F606IF — ES14:W) (yR) ^0.08. An alternative ex¬ 
planation put forth bv IV13l for the offset seen in optical 
colors is the downward revision of the GGG metallicity 
scale, although this would require a shift of [Ee/H] >0.3 
dex based on the relations of Eig. [TTJ This is clearly un¬ 
reasonable given the quality of the spectroscopic metal- 
licities (e.g. Eig. [HI. 

The color offsets seen in the case of the DSED models 
are less drastic, and are confined to absolute values of 
A{J — Ks) which fall well within the margin allowed by 
photometric error and calibrat ion u ncertainties. While 
it may be unsurprising that the lDlOl models require min¬ 
imal color shifts since we adopted distances and redden¬ 
ings based on the fits of these same models at optical 
wavelengths, our results serve as the first test of these 
models in the 2MASS filter system over the range of (age 
and metallicity) parameter space occupied by GGGs. Eor 
both the DSED and VR models, the largest (absolute) 
offset is seen for NGG 5927, the most metal-rich cluster 
in our sample. If we recalculate t he off sets for this clus¬ 
ter assu ming the value of [Ee/H] (jDlQ[ )=-Q.5 rather than 
[Fe/i7] (jGQ9( )=-Q.27 (shown as squares in Eig. [2Q|) . the 
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NGC0104 
[Fe/H]=-0,76 
DSED;-0.001 ±0.011 
VR:-0.012*0.011 


NGC6496 

[Fe/H]=-0.46 


NGC6304 

[Fe/H]=-0.37 


NGC5927 

[Fe/H]=-0.29 



Fig. 20.— Cluster photometry shifted to the absolute plane as in Fig.[7l with radial cuts applied in some cases as in Fig. [5] DSED (blue) 
and VR (red) models are overplotted, but after applying the fixed mean color offset from Table O These mean offsets and their standard 
deviation are given in each panel for each model. 


discrepancy between models and data is brought into the 
range occupied by the remainder of the GGCs, but it is 
not resolved. Although our intention is not to sanction 
individual values for cluster parameters, such a down¬ 
ward revision of the NGG 5927 metallicity is suggested 
by the fits of Figs. [TTl [12] and, to a marginal extent. 
Fig. [ITl Also, Fig. [20] suggests that the DSED and VR 
models have essentially identical main sequence colors at 
fixed metallicity, and would continue to fall redward of 
the observed cluster MS in the near-IR at the extremes of 


the [Fe/H] range sampled, although deeper photometry 
is needed to investigate this effect quantitatively. Lastly, 
it is possible that the color offsets between models and 
data could be resolved by a shift in cluster distances. 
However, throughout our investigation we have charac¬ 
terized the difference between models and data in terms 
of a color offset rather than a magnitude offset due to the 
verticality of the RGB in the near-IR. For example, even 
a modest color offset of A(J — Ks) ^0.02 mag would 
require a shift of ^0.2 mag in the GGG distance scale. 
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[Fe/H] (C09) 

Fig. 21.— Fixed mean color offsets applied to VR (gr ey open cir¬ 
cles) and DSED models (black filled circles) in Fig. 1201 shown as a 
function of cluster [Fe/H] from lCO^ . Vertical error bars represent 
the standard deviation of the color difference between the model 
and observed ficucial sequence from Table [5] Offsets for NGC 5927 
calculated using [Fe/if](D10)=-0.5 rather than [Fe/H](C09)=- 
0.27 are shown as open and filled squares for VR and DSED models 
respectively. The horizontal solid line represents equality, and hor¬ 
izontal dashed and dotted lines represent the median and median 
absolute deviation, respectively, of the color offset across all target 
clusters. 

7. SUMMARY AND CONCLUSIONS 

We have presented 2MASS-calibrated J, Ks CMDs and 
fiducial sequences for 12 GGCs. These fiducial sequences 
have been used to produce relations between photomet¬ 
ric indices which describe the shape of the upper RGB 
versus cluster metallicity, in terms of both [Fe/H] and 
[M/H]. The resulting relations have slopes in excellent 
agreement with previous studies, and show zero point dif¬ 
ferences which can be attributed to uncertainties in the 
distance and reddening of the target clusters. While we 
have chosen to use c luste r distances and reddenings based 
on isochrone fits by IDIQI to deep optical photometry, our 
relations are largely insensitive to current uncertainties 
in cluster distances, reddenings and metallicities for the 
more controversial cases. However, their precision as well 
as the size of the present sample could be improved by de¬ 
tailed spectroscopic abundances for a statistically repre¬ 
sentative quantity of cluster members, especially among 
the most metal-rich GGGs. 

A comparison of empirical fiducial sequences to five 
sets of evolutionary models reveals that DSED and VR 
models most successfully reproduce the observed mor¬ 
phology of the RGB, although a color offset is needed to 
reconcile models and data. Although this offset is not 
formally significant compared to photometric errors as 
well as photometric calibration uncertainties (at least in 
the case of t he DS ED mod els), a comparison between 
the results of IDIQI and IV13l suggests that it is unlikely 
to be caused solely by uncertainties in cluster distance 
and reddening values. Models also suggest that a mod¬ 
est {AY =0.04) helium enhancement negligibly affects the 
RGB morphology, although (^-enhancement plays a sig¬ 
nificant role, particularly on the upper RGB. This is in 
accord with our empirical relations, which generally give 
decreased rms residuals when fits are performed versus 
the global metallicity [M/H] rather than [Fe/H]. 


Relations between cluster metallicity and the near-IR 
magnitudes of the RGB bump and the HB (in the case of 
clusters with sufficiently red HBs) are also insensitive to 
the choice of recent cluster distances, reddenings, metal¬ 
licities and ages beyond current uncertainties. Impor¬ 
tantly, when using the LE peak to characterize the HB 
magnitude and its uncertainty, we find a non-negligible 
dependence between the near-IR HB magnitude and the 
cluster metallicity. This dependence, at the level of at 
least SMk{HB)/S[M/H]=-0.4 mag/dex depending on 
the choice of cluster par ameters, is larger than th at pre¬ 
dicted by the models of iSalaris fc Girar'^ (j2QQ2D and is 
robust to the choice of cluster distances, reddenings and 
ages. 

The photometric catalogs and observed fiducial se¬ 
quences presented here are being made publicly avail¬ 
able, so that the our relations may be modified as the 
GGG distance and metallicity scales are improved. Sim¬ 
ilarly, it is our hope that the 2MASS-calibrated photo¬ 
metric catalogs may be of use as secondary standards 
for near-IR adaptive optics imagers on large telescopes, 
where saturation and/or a small field of view can impede 
photometric calibration. 
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